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Abstract. - We study the effects of crosscorrefations of noises on the scaling properties of 
the correlation functions in a reduced model for Magnetohydrodynamic (MHD) turbulence [A. 
Basu, J.K Bhattacharjee and S. Ramaswamy [Eur. Phys. J B 9, 725 (1999)]. We show that 
in dimension d crosscorrelations with sufficient long wavelength singularity become relevant 
and take the system to the long range noise fixed point. The crosscorrelations also affect the 
ratio of energies of the magnetic and velocity fields (Et,/E v ) in the strong coupling phase. 

^sj ■ In dimension d = 1 the fluctuation-dissipation theorem (FDT) does not hold in presence of 

short range crosscorrelations. We discuss the possible effects of crosscorrelations on the scaling 

yp , properties of fully developed MHD turbulence. 

o 

On 

Numerical studies of MHD turbulence in steady state demonstrate the existence of scaling 
(^) • and multiscaling properties of the structure functions jjj, |2j which are different from fluid 

f (b(k).b(-k)) 

turbulence. The ratio Eb/E v = *£ is believed to be an important parameter in 

J k (v(k).v(-k)) 

characterising statistical properties of MHD turbulence: E\,jE v = =>■ fluid limit, E\,jE v << 
1 => the kinetic regime, Eb/E v ~ 1 =>• the equipartition regime, and Eb/E v >> 1 => the 
magnetic regime (dominated by a mean magnetic field) [Q. As Eb/E v increases from to 1, 
one should be able to observe a crossover from fluid-like behaviour ro MHD-like behaviour - this 
has been confirmed in a recent shell model study |2|,[|. In a model for MHD turbulence Et/E v 
appears as a derived quantity depending upon the viscosities and the external forcings. So for 
a theoretical understanding, it is important to know which parameters influence E\,jE v . MHD 
turbulence is governed by the equations of magnetohydrodynamics (3dMHD) 01. However, 
for simplicity we work with the one-dimensional (Id) reduced model of MHD jp| and its d- 
dimensional generalisation. In this letter wc analytically examine the scaling properties of the 
model equations, driven by stochastic noises. We, in particular, study the effects of noise- 
crosscorrelations. Our main results, obtained by applying renormalisation group and self- 
consistent mode coupling methods on the model equations (see below) include the existence 
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of a roughening transition for d > 2 for short range noises, the dependence of Ef,/ ' E v on the 
crosscorrelations in the rough phase, and breakdown of the FDT in Id. The model equations 
can be easily extended to d-dimension by considering ^ — > V and u, b as vectors (see Ref. H 
for a study of decaying MHD turbulence using the d-dimensional version of the model of 

Ref. |): ' 

f + ^V„* + ^ = ,VW, (1) 

ah 

— + A 3 V(u.b) - M V 2 b + g, (2) 

at 

where v and /x are fluid and magnetic viscosities respectively, f and g are external forcing 
functions. The Galilean invariance of the Eqs.([l]) and (||) requires that Ai = A3; A2 can be left 
arbitrary M . The external forces f , g are taken to be stochastic, with zero mean and Gaussian 
distributions. We choose, 

(/i(k,t)/ i (-k,0)) = 2k i k i D l {k)5{t), (3) 

( 9i (k,t) gj (-k,0)) = 2kikjD 2 (k)6(t), (4) 

C/Kk,i) ft (-k,0)) = 2D t3 (k)S(t). (5) 

Care must be taken while fixing the structure of Dij(k): It should be imaginary and odd in 
k H due to the specific parity properties of u and b. We break up Dij into its symmetric 
and antisymmetric parts (it is easy to see that in Id only the symmetric part survives): 
D^ = DL + D" with £)?■ = ikikjD(k) with D(k) = — D(— k) is the symmetric part and we 
set the anitisymmetric part Df, = for time being. Then the Eqs.p] and @ can be converted 
to those of Erta§ and Kardar, which describe the properties of drifting polymers in a medium 
by the simple transformations u = Vh, b = V(j> and f u = V771, fb = V772 

f + y(VM 2 + ^W) 2 = ^ + , 1 , (6) 

^ + A 3 (V/i).(V</>)= Ai V 2 </> + 7 ?2 . (7) 

We choose D x (fc) = D x + D p k~ p , D 2 (k) = D 2 +D s k- S , with D 1 ,D 2 ,D p and D s being con- 
stants. We also assume |-D(k)| ~ Dk~ y . We are interested in the long-time, long-wavelength 
properties of the model; i.e., we want to calculate the two roughness exponents Xu and Xb of 
u and b and the dynamic exponent z. We employ a standard dynamic renormalisation group 
procedure. Due to Galilean invariance, none of the nonlinearities will renormalise |5(]. Thus, 
if we were to carry out a renormalisation-group transformation by integrating out a shell of 
modes Ae~ £ < q < A, and rescaling r — ► e r, u — ► e e%u u, b — > e lXb h, t — > e iz t, the couplings 
Ai and A2 would be affected only by the rescaling: Ai — » e i ( Xu+z ~ 1 ^\\, \ 2 — » e e ( 2xb ~ Xu+z ~ 1 * > \ 2 , 
implying Xu = Xb = 1 ~ z - These are true independent of the forms of Di(k), D 2 (k) and -D(k). 
Note that the roughness exponents Xi an d X2 of the fields h and (/> are related to those of u 
and b by Xu = Xi - 1, Xb = X2 - 1- 

At the beginning let us assume D p = D s = D = 0, i.e., the noises are short ranged 

and there is no crosscorrelation. We define two dimensionless coupling constants U = A Dl 

and V = A ^ 2 . The nontrivial fixed point is given by U — V — 2d I 3 which is stable 
at d = 1 implying a rough phase for both h and (f> correlations and unstable at d = 2 + e 
indicating a smooth-to-rough transition for both h and <j> correlations; this is exactly the KPZ 
universality class and the fixed point is the standard KPZ fixed point. At the fixed point 



Abhik Basu: Noise-crosscorrelations and Magnetohydrodynamics 




* * 




Fig. 1 - One-loop diagrams contributing to the renormalisation of D. Each of them is individually 
zero (see text). 

v\ = v%, D\ = D%, giving in Id, z = 3/2, Xu = Xb = —1/2 and in dimensions d = 2 + e, 
z = 2 + O(e) 2 , x« = Xb = — 1 — 0(e 2 ). These results imply that at the fixed point the model 
equations decouple in terms of the Elsasser variables z ± — u±b each of which obeys usual 
Burgers equation with (f + .f~) = (f + ,f~ are the noises in the equations for z + ,z~), i.e., 
noises also decouple in this new representation. Next we make D p and D s nonzero. In Id 
p, s < 1/4 the long-range part of noise correlations are irrelevant - this is very similar to 
Ref. [pi (see also Ref. [pi for a detail discussion). However when any of p, s > 9 situation 
changes drastically. It is usually taken that the long-range component of the noise does not 
get renormalised under the application of RG. It is however, easy to see that if p 7^ s this will 
no longer be true; the less singular one will be renormalised to become as singular as the other 
one. At the fixed point we again find v = \x. Using Galilean invariance and nonrenormalisation 
conditions on the long-range components of the noises the exponents can be calculated exactly 
in dimension d ||,g]: We find 

d + 1 - 2p 2p-d-l 
z = H jj , Xi=X2= g ■ (8) 

At this long-range noise fixed point the model again decouples into to two Burgers equations in 
terms of the z variables as defined above. However the noise correlations do not necessarily 
decouple as in general D p ^ D s , i.e., (/+ /~) ^ 0. Let us now examine the scaling properties 
in presence of a nonzero crosscorrelation, i.e., D ^ 0. We first consider the case when D p = 
D s = 0. The procedure of one-loop calculations remain same. There are one-loop diagrams 
(see Fig{y) which in principle renormalise D, however their numerical values are zero, i.e., D 
does not renormalise. Each of the diagram has the structure 

dQd d q D 1 ( q )iD(q) 

2tt 2tt (fi 2 + vq^f ( ' 

as the integral is odd in q (since -D(q) is a pseudo-scalar and odd in q). Due to a non zero D 
there are now additional contributions to both D\ and D2 (see Fig.0), with same magnitude 
but opposite signs, clear whether D\ = Di at the RG fixed point. However, recursion 
relations for v and /1 remain unchanged. Hence at the RG fixed point v = fi. We define a 
new dimensionless coupling constant U = X J? 2 in addition to U and V which we have defined 
earlier. The flow equations for the couplings are 
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Fig. 2 - A figure showing additional one-loop diagrams renormalising (a) D\ and (b) D2 respectively. 



and 



dU 



d + 2y + 3U- 



d-1 



W- 



d-2. 



(12) 



Ad Ad ' 

Let us analyse the case when the crosscorrelations are short ranged, i.e., y = 0. It is easy to 
see that the nontrivial fixed point solutions of the Eqs.10 ll],[l2| are the same as the case when 
there were no crosscorrelations: U = V — 2d -3 an( ^ U — 0, representing again a rough 
phase at d = 1 and smooth-to-rough transitions at d > 2. Interestingly these fixed points are 
marginal in U in a linear stability analysis. All these suggest that a two loop calculation is 
necessary to calculate the full phase diagram. In Id this holds good till y crosses the critical 
value 1/4 as given above. As y becomes larger than 1/4, system crosses over to the long range 
noise fixed point with the exponents given by the exact values (in any d) 



1 



l-2y 



Xu 



Xb 



2y-d-l 



(13) 



Note that at the one-loop level fluctuation corrections to D\ due to the crosscorrelations (i.e., 
D) is positive, i.e., (uu) autocorrelation is enhanced due to D, whereas it is negative for 
Z?2, i-e., (bb) autocorrelation is reduced from its bare value. Since b(r, t) is real (bfeb_fc) is 
positive definite. This suggests that in this model, for positive definiteness of (bfcb_fc) one 
cannot have arbitrarily large D. This property is also shared by 3dMHD equations |14|. If 
one also introduces long range autocorrelations, i.e., one has nonzero D pi D s in the problem 
then when all of y, p, s are less than 9 the KPZ fixed point is stable. However if any of y, p, s 
is larger than 9 then the system crosses over to the long range noise fixed point like before 
(for positivity of (bfcb_fc) we must have p > s in this model. If all or more than one of 
them is larger than 9 then the exponents are goverened the largest of them. Here also an 
inetersting point is that if y is larger than p or s then both D p and D s renormalise such that 
renormalised auto correlations scale as k~ y ; however if y < p, s then crosscorrclation does not 
renormalise (this is again due to the oddness of the crosscorrelation). All these suggest that 
crosscorrelation can be a relevant operator in a a field theory for this problem. 

Having considered the weak coupling limit for d > 2 we now examine the strong coupling 
phase by using a one-loop selfconsistent scheme. In particular we investigate the depen- 
dence of the ratio D2/D2(^ Eb/E v ) on the crosscorrelations. For simplicity we work with 
short-ranged noises only. We also calculate the upper critical dimension (d c ) of the model 
in presence of short range crosscorrelations. The model is only logarithmically rough above 
d c . The value of d c is still is not well settled for the KPZ equation. Some workers sug- 
gest it to be 00 p| ; some others suggest a finite d c : Lassig and Kinzelbach |12[ mapped 
KPZ equation onto a problem of a directed polymer in a random medium and showed that 
d c < 4; Bhattacharjee p3[ used a mode coupling approach to show that d c — A. We fol- 
low Ref. E3| closely to find whether d c depends upon D. For simplicity we work with the 
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height variable representation of our model equations, i.e., with Eqs.o and u[ We assume 
scaling forms G h = k~ z g h {oj/k z ); C h = k~ d - 2 ^- z f h {u/k z ); G = k' z g^uj/k z ); C = 
k -d-2 X 2-zf g ( u /k*); C = Z)(k)fc- d -xi-X2-^ f(u)/k z ), with D(k) = --D(-k) and |Z)(k)| = L>, 
a constant, i.e., C is a pseudoscalar in k. Due to Galilean invariance of the nonliearities 
Xi — X2 = X- From diagramatics, we have seen that one-loop corrections to v and \x are 
same. Hence we put v — /i without any loss of generality. This implies Si(k, uj) = £ 2 (k, uj) = 
S(k, uj). We assume the forms for the zero-frequency response and correlation functions as 

£(k,w = 0) = Tk z , Ci(k,w = 0) = L>ifc- 2x - d - 2 , C 2 (k,w = 0) = D 2 k' 2x - d - z . (14) 

We employ a small x expansion as used in Ref . [Jl3| . We calculate self-consistent expressions 
for the self energies and correlation functions: Matching at uj = pi 

_ll_ _ I f A [(I-P)-Pl 2 f1 , Dl , ^ 2 
Z)^ 2 2 J {2n) d p d + 2 *(k - p) d + 2 x { Z> 2 D 2J ' 

r 2 y d d P [(i-p)-p] 2 ,A #\ n ^ 



L> 2 A 2 J (27r) d p d + 2 x(k-p) d + 2 x v D 2 £> 2 

Since there is no one-loop diagrammatic correction to C there is no self energy correction for 
C. To the leading order, from one-loop self energy one obtains 

r 2 Sd i, 1JL ^»2s r 2 s d i n ,D! 



DiA 2 {2ir) d 2d y D t " D 2 X 2 {2n) d 2d y D 2 ' 

where Sd is the surface of a <i-dimensional sphere. On the other hand, one-loop expression 
for correlation function gives, by extracting the high momentum parts (p >> 1) following 
Ref. jT§ 



DiA 2 4 (27r) d p d - 2 +3x L ^ y y D JI, D 2 X 2 4 (2tt)<V- 2 + 3 x l L> 2 v L> 2 
From EqsJlT] we find 

^1 + 2^/3 + 2/3-1 = 0, (18) 

•i-'l -L'2 

where /3 = (-jr) 2 - Notice that for = jy- = 1. Since one-loop Eq.18 is correct upto 0(/3), we 



look for solution for jf- of the form -^ = 1 + a/3, such that for /3 = we recover £>! = D 2 . We 
obtain a = —2, i.e., D 2 /D\ = 1 — 2/3. So within this approximate calculation /3 cannot exceed 
1/2 (i.e., D < Di/2 ). This immediately gives Xi = X2 = 1/2 + 0((i 2 ) in ci = 1 in agreement 
with the one loop DRG results: xi = X2 = 1/2 which are correct to 0(0). Following Ref. [jl3| 
we obtain d c = 4 + 0(/3) 2 . These results also suggest the need of a two-loop calculation. 

We now discuss the effects of the antisymmetric part of the crosscorrelation (we now work 
with Eqs.[HE]as, in the presence of antisymmetric crossocrrelations, Eqs.[t],|| do not reduce to 
Eqs.|| and|7]. However, the procedure is identical.). We have 

D£(k) = -D£(-k) = D£(-k) = -[D£(k)]*. (19) 

We choose Dfj(k)D° , i (—'k) = D 2 k 4 and D = 0. When the crosscorrelation is antisymmetric 
(i.e., D = 0, D + 1 0) fluctuation corrections to both D\ and D 2 are positive and same value: 
The flow equations for dimensionless coupling constants U, V and U = ^-p- are 

dU Tr/n n rr2 2d-3 V 2 U 2 
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and 

f-»p-« + * + »^ + «'^]. (22, 

Since the fluctuation corrections to Z?i and Z?2 due to D are positive, there is no restriction 
on the strength of D, unlike the case with symmetric crosscorrelations. The exponents at 
the weak coupling fixed point for d > 2 remain unaffected by U at the one loop level (at 
Id, U = 0). It is easy to find out d c when D ^ following Ref. @. We find D x = D 2 at 
the strong coupling fixed point and d c — 4 + 2/3, keeping terms only at the linear order in /3. 
Thus in presence of an additional noise in the form of the crosscorrelations, d c and hence the 
roughness (at d c the roughness exponent drops to zero from a positive value) increases. 

We have seen that when D = (i.e., there are no crosscorrelations) D\ = D 2 (or U = V, 
by assuming v = /i) at the RG fixed point. In fact, U — V is the FDT line in the (U — V) 
plane. If we start with D\ = D 2 = D, this remains invariant under the RG transformations. 
Along this line the model equations decouple into two independent Burgers equations in terms 
of z + ,z~. In Id along the line U = V vi/v = Di/D, telling immediately that the Fluctuation- 
Dissipation Theorem (FDT) pd holds. Thus we get the steady state probability distribution 
P[u, b] ~ exp{jj f x {(d x it) 2 + (d x b) 2 }]. We immediately obtain Xu = Xb = —1/2 exactly and 
hence z = 3/2. However in Id when D ^ the exponents do not depend upon D at the 
one-loop level, FDT does not hold any longer, simply because there are additional fluctuation 
corrections to D\ and D 2 due to D: 



A 2 L>i \ 2 D 2 , , \ l D Y \ 2 D 2 , , , 

* = "(1 + -^ + -^), /» = Md + -^ + ^), (23) 

^=A[l + ^ + ^^ + 2^(#) 2 ],^=^l + 2^-2^ ( A )2 a )2] . 

1 4^ 3 4L>i M 3 4^ 3 V Di J ' " L 4^ 3 4^ 3 V Di V L>i J 

(24) 
Obviously vijDu ^ v Q jD\ and D%/Di = 1 — 2/3^0, consequently the FDT does not hold 
good. Thus £>i = Z?2 is maintained under RG transformations. Hence the exact values of 
z j Xi) X2 are not known, unlike the case with D = 0. Recall that when £> = 0, vi/Du = v/D\, 
which is nothing but a statement of the FDT [^O). Note that the dimensionless ratio {jy) 2 
determines the amount of deviation from the FDT. 

What are the implications of the above mentioned results on fully developed MHD turbu- 
lence. In MHD turbulence in general crosshelicity is nonzero and found to scale as ~ fc" 5 / 3 
(within error bars) fcj. This makes it imperative to have nonzero crosscorrelations of noises 
in stochastically driven MHD models (only the symmetric part contributes to the crosshe- 
licity, however, the existence of an antisymmetric part cannot be ruled out on the basis of 
any general principle). fc~ 5 ' 3 -spectrum requires noise correlation to scale as fc~ 3 in 3d (for 
both auto and crosscorrelations). If all bare correlations scale in the same way then one loop 
corrections to the autocorrelations at zero external frequency scale as the bare ones in the 
low momentum limit; however one-loop correction to the crosscorrelation is identically zero. 
This means the energy spectrum is at least as singular as (or more than) the crosshelicity 
spectrum. It is interesting that similar behaviour holds good also for the stochastically driven 
3dMHD equations E3]. We have also seen that in our calculations, the ratio D 2 /Di depends 
upon D, D, the strength of crosscorrelations. Now E-^jE v ~ D 2 /Di ~ 1 — 2/3. Thus in this 
model, by varying D and D we can achieve different values of Eb/E v , ranging from to 1. 
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Similar dependence of Eb/E v on D holds for 3dMHD equations also E3. Of course, to model 
MHD turbulence one must work with noise correlations that are singular in the k => limit. 
For such noises, scaling (i.e., the fc-dependence) of the noise correlations do not change under 
renormalisation, but the amplitudes change in such a way that D\ is enhanced and Di is 
suppressed in presence of a crosscorrelation, affecting the ratio Eb/E v . Thus, in conclusion, 
we have discussed the role of a crosscorrelation of noises on the scaling properties of the re- 
duced model for MHD. In particular we have shown that sufficiently singular crosscorrelation 
can be a relevant perturbation on the model. The upper critical dimension of the model in 
presence of short range noise only does depend upon the strength of the crosscorrelations and, 
interestingly, the dependence is different for the symmetric and the antisymmetric parts. In 
the presence of antisymmetric crosscorrelations, the reduced model does not reduce to that 
of Ertas and Kardar. However one cannot rule out its presence in the context of MHD on 
any symmetry grounds. It is important to examine the effects of crosscorrelations (thus finite 
crosshelicity) on the values of the universal numbers in fully developed MHD turbulence, e.g., 
Kolmogorov constants, intermittency exponents |[L4| . It would be very interesting to check 
our results in numerical simulations of our model and/or 3efiVIHD equations. 
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